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ABSTRACT 

Q i We present Gemini-North GMOS/IFU observations of a young star cluster and its 

' environment near the centre of the dwarf irregular starburst galaxy NGC 1569. This 

C/3 , forms part of a larger and on-going study of the formation and collimation mechanisms 

of galactic winds, including three additional IFU pointings in NGC 1569 covering the 
I base of the galactic wind which are analysed in a companion paper. The good spatial- 

,—1 . and spectral-resolution of these GMOS/IFU observations, covering 4740-6860 A, allow 

^ ' us to probe the interactions between clusters and their environments on small scales. 

0^ ! For cluster 10, we combine the CMOS spectrum with HST imaging to derive its 

■ properties. We find that it is composed of two very close components with ages of 5- 
' 7 Myr and ^5 Myr, and a combined mass of 7± 5 x 10'^ M0. A strong red Wolf-Rayet 

(WR) emission feature confirms our young derived cluster ages. 
00 ' A detailed analysis of the Ha emission line profile shapes across the whole field- 

' of-view shows them to be composed of a bright narrow feature (intrinsic FWHM 

, ~ 50 kms~^) superimposed on a fainter broad component (FWHM ^ 300 kms~^). 

■ By mapping the properties of each individual component, we investigate the small- 
^ I scale structure and properties of the ionized ISM, including reddening, excitation 

•rH , and electron densities, and for the first time find spatial correlations between the 

' line component properties. We discuss in detail the possible mechanisms that could 

, give rise to the two components and these correlations, and conclude that the most 

■ - - i likely explanation for the broad emission is that it is produced in a turbulent mixing 

layer on the surface of the cool gas clumps embedded within the hot, fast-flowing 
cluster winds. We discuss implications for the mass-loading of the flow under these 
circumstances. The average radial velocity difference between the narrow and broad 
components is small compared to the line widths, implying that within the IFU field- 
of-view, turbulent motions dominate over large-scale bulk motions. We are therefore 
sampling well within the outer bounding shocks of the expanding superbubbles and 
within the outflow 'energy injection zone'. 
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1 INTRODUCTION 

This is the first in a series of papers presenting integral field 
unit (IFU) observations of ionized gas in the dwarf star- 
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burst galaxy NGC 1569. The study is aimed at examining 
the small-scale relationship between the interstellar medium 
(ISM) and star-cluster population. In this paper we focus on 
the data reduction techniques and present the first results. 

At a distance of 2.2 ± 0.6 Mpc (|lsraelll 19881 ). NGC 1569 
is one of the closest examples of a starburst and an excellent 
analogue to high-redshift dwarf starbursts implicated in the- 
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ories of galaxy for mation. The total B magnitude of NGC 
1569, Mb = -17 jlsraeJIlQSSl}^ and a remarkably homoge- 
neous metallicity of 0.25 Zp, (jPevost. Roy fc Drissenlll997l . 
iKobulnickv fc Skillmanlll997l ) makes this galaxy similar (in 
some ways) to the Small Magellanic Cloud. 

NGC 1569 has attracted considerable attention because 
observations at a variety of wavelengths show that it has re- 
cently undergone a galaxy-wide burst of star formation that 
peaked between 100 Myr and 5-10 Myr ago (|Grege:io et al.l 
1 19981 ). The presence of numerous Hll regions reve als that 
there is still substantial ongoing star formation (IWalleJ 
Il99ll). T he two well-known super star clusters (SSCs) A and 
B (|AbIes|[l97ll : lArp fc Sandagelll985l ') are prominent prod- 
ucts of the starburst episode, which may have been triggered 
by an interaction w ith a nearby H I cloud (|Stil fc Israelii 199^ : 
iMiihle et al.ll200"5l ). An extended syste m of Ha filaments, in- 
dicative of a bipolar outflow is seen (jHeckman et al.lll995l : 
iMartinI Il998l ). X-ray observations show that the wind is 
metal-enr iched and emanates from the full ext ent of the 
Ha disc l|Martin. Kobulnickv fc Heckmaril |2002[). but only 
from ~0.6 times the extent of the H i disc (|Walle"il Il99ll : 
IStil fc Israe3l2002l : iMartin et al.ll2002l . their fig. 7). 

The star-cluster population of NGC 1 569 has been 
studied extensively using HST observations. iHunter et al.l 
(|2000l ) identified and catalogued a total of 48 compact but 
resolved clusters using W FPC2 imaging . Usin g the same 
data but different criteria, lAnders et al.l (|2004l 'l found 179 
clusters, with the majority being formed in an intense 
starburst event which began ~25 Myr ago. Of the clus- 
ter population, only the mo st luminous cluster s A, B and 
No. 30 (nomenclature from iHunter et al.l |2000| ) have been 
studied in any detail. These three SSCs alone provide a 
significant fraction (20-25 per cent) of the total optical 
and near-infrared l ight in the central region of NGC 1569 
jOrigl ia et al.ll200l|). Cluster A has actu ally been resolved 



J lex, ionized structure with filame n ts, bubbles and loops 
Hunter. Hawlev fc Gallagheil Il993l: iTomita. Ohta fc Saitd 



into two clusters ( de Marchi et al ] |l997h . and has an in- 
tegrated age of ~7 Myr, with a prob able small age dif- 



ference between t he two components (|Hunter et aLl lioOOl : 
iMaoz et all 1200 ll : lOiTglia et al.l 120011 ). Clusters B and 30 
are older with ages of 10-20 and ~30-50 Myr respectivel y 
jHunter et al.ll2000l : lOriglia et allbOOll : lAnders et al]|2004 ) . 

The most conspicuous active star-forming region is lo- 
cated 90 pc to the west of cluster A and c orresponds t o 



the brightest Hll region in NGC 1569 (No. 2: IWallej|l99 j) 
and the peak of the ther mal radio emission (|Lisenfeld et al.l 
I2OO4I ). the dust emi ssion (iLisenfeld et al ll20oi ) and the sub- 
mm/FIR emission (|Galliano et al.ll2003l ). This region is at 



the e astern edge of a large CO cloud complex (|Tavlor et al.l 
1 19991 ). which contains an d is surrounded b y a number of 
radio continuum sources (|Greve et al.ll2002l ). IHunter et al.l 
(|2000l ) find two clusters in this region, numb e rs 6 a nd 10, 
the latter of which is found by iTokura et al.l (|2006h to be 
surrounded on one side by mid-IR [S iv]10.5 ii.m, emission 
coincident with a gas knot clearly visible in Ha. Cluster 10, 
one of the subjects of this paper, is the third v isually bright- 
est cluster after A and B (|Hunter et al.|[200ol ). but has not 
been studied in any detail thus far. 

Stellar-wind and supernova-driven outflows powered 
by the collective injection of kinetic energy and momen- 
tum from massive stars in starbursts can drastically af- 
fect the structure and subsequent evolution of dwarf galax- 
ies. In NGC 1569, Ha images show a chaotic, com- 



I1994I : iHeckman et al.1 1 19951 : lAnders et al.1 12004) w h ile th e 
neutral ISM is highly turbulent (jStfl fc Israeli I2OO2I ). 



iHeckman et al.l (|l995 ) observe Ha line profiles which are 
"distinctly non-Gaussian", exhibiting broad, asymmetric 
wings that they decompose using multi-Gaussian profile fit- 
ting. Within 2" of SSC A, the Ha line has a full-width 
zero-intensity of 30-50 A and they find the wings to ac- 
count for up to ~30 per cent of the total Ha fiux within 
this region. Broad emission line wings have been detected 
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Westmoquette et al. (2007b, in prep.) present high- 
resolution _ffST/STIS observations of the central regions of 
M82, and identify a ubiquitous underlying broad component 
to the optical nebular emission lines in this galaxy. They 
were able to track the properties of this component through 
the central starburst regions, and relate them to the highly 
fragmented state of the ISM. 

The nature of the energy source for these broad lines is 
not clear, and a number of contesting explanations have been 
proposed, all of which relate to the action of stellar winds 
and supernovae (SNe): broad lines resulting from the hot, 
turbulent gas confined to large bubbles blown by the winds; 
large velocities associated with supernova remnants (SNRs) ; 
integrating over many ionized structures (shells) with differ- 
ent expansion velocities; or from blow-out phenomena such 
as champagne flows or superwinds. These mechanisms act 
to differing degrees over differing spatial extents depend- 
ing on the ambient physical conditions; the line widths and 
integrated shapes thus change accordingly. This makes an 
assessment of their effect or importance difficult, particu- 
larly since the available observations are often conflicting. 
Through the data presented here, we attempt to clarify 
which mechanisms could apply in the case of NGC 1569. 

HST narrow-band Ha images, in particular, illustrate 
that the ISM contains many small-scale structures (e.g. 
iBuckalew fc Kobulnickvl l2006h . To probe the roots of the 
wind outflow from NGC 1569, it is necessary to investigate 
the interaction of the individual winds from clusters with 
their environments on high angular-scales. Good spectral- 
resolution and spatial coverage are, however, equally impor- 
tant for determining properties such as gas dynamics and 
excitation mechanisms. Integral fleld spectroscopy (IPS) is 
ideally suited to these type of observations, and modern, 
large format integral field units (IFUs), such as those found 
on 8m-class telescopes, currently provide the best opportu- 
nity to fulfill these demanding requirements. Historically, it 
has been very difficult to manage the data products from 
such instruments, but in recent years the infrastructure to 
deal with the reduction, analysis and visualisation of IPS 
data has improved dramatically. 

To this end we have obtained Gemini-North 
GMOS/IFU observations of the cluster wind-ISM in- 
teraction zone in the central region of NGC 1569. Our 
observations focus on the small-scale and are of high 
enough spatial-resolution to study the wind-ISM inter- 



A young star duster's environment in NGC 1569 3 



action processes in detail. In this first paper we present 
the dataset for one of four IFU positions, and describe 
in detail the reduction and analysis techniques employed 
and the software tools we have developed. In a second 
paper (Westmoquette et al. 2007a; Paper II) we present the 
analysis of the remaining IFU positions covering the central 
outflow region. In a third paper (Westmoquette et al., in 
prep.; Paper III) we will examine the properties of the outer 
wind regions using deep Ha imaging and IFU observations. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

In November 2004 queue-mode observations using the 
Gemini-North M uhi-Object Spectrograph (C MOS) Integral 
Field Unit (IFU: lAllington-Smith et al.ll2002l ') were obtained 
covering four regions near the centre of NGC 1569 (pro- 
gramme ID: GN-2004B-Q-33, PI: L.J. Smith), with 0.5- 
0.8 arcsec seeing (6-9 pc at the distance of NGC 1569). 
A nearby bright star was used to provide guiding and tip- 
tilt corrections using the GMOS on-instrument wave front 
sensor (OIWFS). 

Depending on the combination of spatial and spectral 
coverage required, the GMOS IFU can be operated in one- 
or two-slit modes. For our purposes, we opted for the one- 
slit mode giving a field-of-view of 5 x 3.5 arcsecs (which 
corresponds to approximately 50 x 35 pc at the distance of 
NGC 1569) sampled by 500 hexagonal contiguous fibres of 
0'.'2 diameter. An additional block of 250 fibres (covering 
2.5 X 1.7 arcsecs) are offset by 1' from the object field pro- 
viding a dedicated sky view. The coordinates of the IFU 
field centre analysed in this paper were a — 04** 30™ 47! 27; 
5 = -1-64° 51' 03'.'2 (J2000) and the total exposure time, split 
into four separate integrations, was 6720 s. Using the R831 
grating to give a spectral coverage of 4740-6860 A and a 
dispersion of 0.34 Apix"^, we were able to cover the impor- 
tant nebular diagnostic lines of Ha, H/3, [O lll]AA4959, 5007, 
[Nii]AA6548,6583 and [Sii]AA6716, 6731. 

The GMOS spectrograph is fed by optical fibres from 
the IFU which reformats the arrangement of the spectra 
for imaging by the detector. This is composed of three 
2048 X 4068 CCD chips separated by small gaps. In order 
to remove the pixel-to-pixel sensitivity differences, and en- 
able the wavelength and flux calibration of the data, a num- 
ber of bias frames, flat-fields, twilight flats, arc calibration 
frames, and observations of the photometric standard star 
G191-B2B, were taken contemporaneously with the science 
fields. 

The IFU positions were chosen to cover selected areas of 
the disturbed ionized interstellar medium in the inner region 
of NGC 1569, near the roots of the wind outfiow. In this pa- 
per we present data for one position which samples cluster 
10 and its immediate environment. In Fig. [T] we show the 
position of the IFU field on an archive HST F656N image 
(GO 6423, P.I. D. Hunter; see Table [T| obtained with the 
Wide Field Planetary Camera 2 (WFPC2). Cluster 10 is lo- 
cated on the east ernmost part of the ionized complex centred 
on Hll region 2 (|Wallerlll99ll ) and is ideal for the study of 
the cluster-ISM interaction. Just off the south-west edge of 
the IFU field are two tentative 1.44MHz detections M-c and 



M-d of lCreve et al] (|2002l ). To the north of the field, ~25 pc 
from cluster 10, is VLA- 1 6, a n on-thermal radio source also 
reported bv lCreve et al.l (|2002l ') who consider it may be an 
extended low surface brightness supernova remnant (SNR). 

2.2 Data Reduction 

Because of the complexity of the data products from IFU 
instruments, a considerable amount of time was invested 
in learning and developing reduction and analysis tools to 
enable us to extract the most information we could from this 
extensive dataset. 

The field-to-slit mapping for the GMOS IFU reformats 
the layout of the fibres on the sky to one long row containing 
five blocks of 100 object fibres interspersed by five blocks of 
50 sky fibres. In this way, all 750 spectra can be recorded on 
the detector at the same time. 

Basic reduction was performed using the Gemini 
pipeline reduction package (implemented in iRAlf^). There 
are a number of tasks within this package that are designed 
to perform certain functions, and each contain switches to 
allow them to be applied at different stages of the reduction 
procedure. The first stage was to run the GFREDUCE task on 
the flat-fields, twilight flats and arc calibration exposures in 
order to prepare the raw files for reduction, and subtract 
the overscan level and bias frames. GFEXTRACT traces and 
extracts the spectrum of each fibre (using a ±2.5 pixel aper- 
ture in the spatial direction), and by running this on the 
fiat-field image, produces a reliable trace of the position of 
each spectrum on the CCD. We then ran GFRESPONSE on the 
GFREDUCE'd twilight fiats to create a throughput correction 
function for each fibre (an essential step to correct differ- 
ences in fibre-to-fibre transmission). The final step before 
reducing the science data was to extract the 750 arc calibra- 
tion spectra from the arc lamp exposures (using GFREDUCe) 
and establish and apply a wavelength calibration solution 

using GSWAVELENGTH and GFTRANSFORM. 

Running GFREDUCE on the science observations using 
the spectrum trace created from the flat-fleld, the through- 
put correction determined from the twilight flat, and the 
wavelength calibration from the arc exposures, creates a 
data file containing 750 reduced spectra, each one pixel in 
width and ordered by the position they were recorded on 
the CCD (which we will call the 'fibre order'). After extrac- 
tion of the spectra, the x and y spatial unit is now termed 
'spaxel' to differentiate from 'pixel' which refers to the CCD. 

Cosmic-rays were cleaned from the data at this stage 
with the Laplacian cosm ic-ray identification routine LACOS- 
MIC (|van DokkumlbOOll l. Fits to the flux of the [Ol]A5577 
telluric emission line, assumed to remain spatially stable in 
intensity for each observation set, for each of the 750 object 
spectra, showed a significant (--^10 per cent) dip from the 
centre to the edges when plotted against fibre order. We de- 
cided that a second throughput correction was needed to re- 
move this variation. This correction was measured by fitting 

^ The Image Reduction and Analysis Facility (iRAp) is dis- 
tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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Figure 1. i?5T/WFPC2 Planetary Camera image taken tliroug li tlie F656N filter o f the central region of NGC 156 9 showing the pos ition 
of the IFU field. A number of the most prominent star clusters jHunter et al.|[200Cll) and radio continuum sources jGreve et al.ll20o3) are 
labelled in black and white respectively. 



a fourth order polynomial to the intensity curve, inverting 
the profile, normalising it to unity, then dividing this into 
the data using the iraf task imarith. 

Fits to the full-width half-maximum (FWHM) of the 
[O l]A5577 line also showed unexpected variations with fibre 
order. We found that the groups of 50 sky fibres at either end 
of the CCD showed elevated widths compared to the central 
groups. Thus to be safe we extracted only sky fibres in the 
central blocks (i.e. fibres 200-250, 350-400 and 500-550). 
After extracting and averaging these fibres for each science 
frame, the resulting spectrum was subtracted from the rest 
of the spectra using imarith. Because of diffuse emission 
from the galaxy extending past the 1 arcmin separation of 
the object and sky fields, there is the possibility that neb- 
ular emission could contaminate the sky spectra. To check 
the level of such contamination, we compared the flux of 
the Ha line in representative spectra with and without sky 
subtraction. The difference in the flux for each case was 1-2 
per cent and therefore deemed insigniflcant. 



2.3 Flux calibration 

As part of the calibr ation dataset, the flux standard star 
G191-B2B (|Okell 19901 ) was observed using the same instru- 
ment set-up as described in Section [2.21 Reduction of these 



data followed the same procedure as described above for the 
science frames, except that we used the throughput response 
file and wavelength solution determined from the reduction 
of one of the other IFU fields. In order to obtain the opti- 
mum S/N spectrum for the standard star, we summed all 
object spax;els across the field of view using GFAPSUM to 
create one averaged spectrum of the standard star. The sen- 
sitivity function was computed using the GSSTANDARD task 
(which includes a second-order extinction correction relevant 
to the Gemini-N observatory), and the fiux calibration for 
each science frame was performed by the GSCALIBRATE task. 

Final combination of the individual exposures for each 
position was done using imcombine. The resulting data file 
was then formed of 750 reduced, sky-subtracted and flux- 
calibrated spectra, one corresponding to each spaxel in the 
field-of-view (including the sky field). Separation of the sky 
fibres from the data files was achieved using the E3D Vi- 
sualisation Tool jSAnchcz, 2003 ) . the development of which 
was a primary goal of the European Commission's EuroSD 
Research Training Network. Each science file now only con- 
tained the 500 object spectra. 

In order to determine an accurate measurement of the 
instrumental contribution to the line broadening, we selected 
spectral lines from a wavelength calibrated arc exposure that 
were close to the Hq and [Olll]A5007 lines in wavelength. 
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and sufficiently isolated to avoid blends. We then fitted sin- 
gle Gaussians to these lines for all 750 apertures and took the 
average. The instrumental broadening (velocity resolution) 
of the final dataset varies between FWHM = 74 ± 5 kms^^ 
at the blue end (with an average S/N of 4.6 in the contin- 
uum), to 59 ± 2 kms^^ at the red end (with an average S/N 
of 14.5). 



2.4 Differential Atmospheric Refraction (DAR) 
Correction 

When observing the spectrum of an object, its light is re- 
fracted by the Earth's atmosphere by varying amounts at 
different wavelengths. This effect, termed differential atmo- 
spheric refraction (DAR), is therefore a function of wave- 
length (horizontal) and airmass (vertical). An advantage of 
IPS over traditional long-slit methods is that it is possible to 
determine and corr ect the effects of DA R using an a posteri- 
ori procedure (e.g. lArribas et al]|l999l l. In order to correct 
our data for DAR, we first had to convert each dataset to 
standard cube format. This involves changing the way in 
which the data are stored, from one spectrum (wavelength 
vs. ffux) per spaxel (spatial coordinate) to a 2-D ffux image 
at each wavelength point. For the DAR correction proce- 
dure, it is required that the data are sampled evenly in x 
and y, and since the format of the GMOS fibres is hexag- 
onal we had to apply an interpolation routine to resample 
the data into contiguous squares with an equivalent 'diame- 
ter'. A diameter of O'.'IS gave the closest match to the num- 
ber of original spaxels when interpolating from hexagons to 
squares. 

Measurement of the DAR shift can be achieved eas- 
ily and automatically if there is an unresolved point source 
in the field-of-view: in our case cluster 10. With the data in 
cube format (2 spatial axes and one of wavelength) and using 
scripts provided by L. Christensen (ESO), the spatial shifts 
were traced by cutting the cubes into a series of monochro- 
matic images, and fitting a two-dimensional Gaussian func- 
tion to the unresolved star cluster to determine the centroid 
position. A third order polynomial was used to fit this trace 
(with wavelength) through the datacubes, giving the relative 
shifts in RA and Dec for each slice with respect to the first 
(reference) wavelength. The measured offsets were then used 
to shift each wavelength slice by then required amount, and 
were the combined back into a datacube with the original 
grid sampling using a modified version of the iraf STSDA^ 
DRIZZLE task (originally written as an implementation of 
the image combina tion method known as v ariable-pixel lin- 
ear reconstruction: [Pruchter fc Hoolj 120021 '). The result is a 
datacube where every fiux at each point in the object corre- 
sponds to the same position in the image for each wavelength 
point. 

The final stage was to crop the edges of each field in 
order to remove spaxels only containing flux for part of the 
wavelength range. This is an inevitable side-effect of per- 
forming a DAR correction. 



^ STSDAS is the Space Telescope Science Data Analysis System; 
its tasks are complementary to those in IRAF. 



2.5 HST images 

To supplement our studies, we obtained HST ACS/HRC 
and WFPC2 broad- and narrow-band images of the central 
regions of NGC 1569 from the HST archive. These images 
are listed in Table [1] (together with the corresponding pro- 
gram IDs and P.I.) and were calibrated using the "on-the- 
ffy" pipeline system. Dithered images were first registered 
using their WCS (World Coordinate System) information, 
then combined using the iraf imcombine routine in order 
to remove cosmic-rays. 

Pig. [2]introduces the concept of imaging with IFU data. 
On the left is the HST ACS/HRC P555W image of the re- 
gion covered in our IFU pointing. In the centre we show a 
'virtual' F-band image reconstructed from our IFU data by 
convolving the spectrum from each spaxel with a l^-band 
filter function. The similarity between the two images is 
obvious once the resolution of each is taken into account 
(ACS/HRC = 0'.'026 pixeP^ GMOS IFU = O'.'IS spaxeP^ 
but seeing-limited at ~0'.'8). One must keep in mind the 
fact that the ACS F555W filter contains contamination from 
nebular [O ill] emission, which is why a certain amount of 
nebulosity is seen surrounding cluster 10. The right-hand 
plot is of the same field but restricted in wavelength to cover 
the continuum only (in the range 6645-6660 A). Overlaid is 
an outline of the 24 spaxels summed to create a spectrum 
of cluster 10 (Section [31). 



3 CLUSTER 10 

Cl uster 10 is the third visually brightest cluster identified 
by iHunter et al.l l|2000l ) but, as can be seen from Fig. [2] 
(left panel), it is in fact a double cluster. Because of the 
close proximity of the two sources, previous studies (e.g. 
iHunter et ai] bOOOD have treated cluster 10 as one object. 
Fitting the light profiles across the length and breadth of the 
IFU field with a Gaussian profile yields a average FWHM 
of 4.7 spaxels or 0'.'84. This corresponds to the maximum 
seeing disc size at the time of observation, meaning that the 
two sub-clusters are not resolved in our IFU data. 

From the F555W ACS/HRC image, we measure the sep- 
aration of the two cluster 10 components as 0'.'35 (3.7 pc). 
The coordinates of the two sources are given in Table[3l Mea- 
surement of the size of the two clusters was made on the 
ACS/HRC F5 55W image, and was achieved using i shape 
jLarsenl Il999l ') together with the tinytim package (jKristl 
l2004h to correct for the point-spread function (PSF). Using 
a circular Moffat function with a power index of 1.5, and a fit 
radius of 0'.'2 gives an effective radius, ReB = 1.18 ± 0.05 pc 
for the north-eastern component (hereafter cluster lOA; see 
Fig. [2] left panel), and R^a = 0.88 ± 0.15 pc for the south- 
western component (hereafter cluster lOB). Due to this very 
small separation, and the crowding of the field around clus- 
ter 10 with other fainter sources, standard aperture photom- 
etry is not possible (and we cannot measure absolute magni- 
tudes). However if all that are needed is relative magnitudes, 
then small-radius aperture photometry is sufficient. Using a 
radius of 4 pixels for the HRC images (F333W, F555W and 
F814W) and 2.5 pixels for the WFPC2 image (F439W) and 
suitable background annuli, we measure (F555W-F814W) = 
1.12 and (F330W-F439W) = -0.27 for the cluster lOA, and 
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Table 1. HST archive images of NGC 1569. 



Filter 


Camera/CCD 


Plate Scale 


Exposure 


Programme 


P.I. 






(" Pix-i) 


Time (s) 


I.D. 




F330W 


ACS/HRC 


0.026 


220 X 3 


9300 


H. Ford 


F439W 


WFPC2/PC 


0.045 


700 X 16 


6111 


C. Leitherer 


F555W 


ACS/HRC 


0.026 


130 X 3 


9300 


H. Ford 


F656N 


WFPC2/PC 


0.045 


800 X 2 


6423 


D. Hunter 


F814W 


ACS/HRC 


0.026 


130 X 3 


9300 


H. Ford 




0.36 0.72 1.08 1,44 1,8 2.16 Z.52 Z.aS 3.24 




Figure 2. Left: HST / ACS HRC F555W image of the region covered by the IFU pointing. The two components of cluster 10 are labelled. 
Centre: IFU data convolved with a y-band filter function. Flight: flux map in the continuum only (6645-6660 A), showing the position 
of cluster 10 and the spaxels extracted to form the cluster spectrum. North is up and east is left. 



(F555W-F814W) = 0.42 and (F330W~F439W) = -0.90 for 
cluster lOB. We note that since the cluster radii are so sim- 
ilar, using the same aperture for both does not introduce 
significant errors. 

Fig.[3]shows the evolutionary path of a Z = 0.004 Star- 
BURST99 (sb99; iLeitherer et al.| [l999l ') instantaneous burst 
evolutionary-synthesis model with a standard Kroupa ini- 
tial mass function (IMF) in U — B, V — I colour space, to- 
gether with the colours that have been determined for both 
cluster 10 components. Also plotted is the colour of the inte- 
grated flux of the two sources (measure d using an aperture 
equivalent in to the cluster size found bv lHunter et al.ll2000l ) 
to compare to parameters derived from the spectrum, since 
the spectrum is derived from a convolution of the fight from 
both components. SB 99 does not give outputs in terms of 
HST filters, so the observed photometry was converted from 
HST to Johnson magnitudes u sing the conversion factors 
given by iHoltzman et al] (|l995h for the WFPC2 data, and 
ISirianni et al.l (l2005l l for the ACS data. From the diagram 
we can now derive an approximate age and an estimate of 
the reddening for each of the cluster components simulta- 
neously. Assuming a extinction value sufficient to deredden 
each point until it intercepts the model track has the follow- 
ing results. For cluster lOA, a.n Ay of «0.8 [E{B-V) « 0.26] 
corresponds to an age of 50-100 Myr, whereas Ay ~ 1.2-2 
[E{B-V) ^ 0.65-0.40] gives an age of 5-20 Myr. For cluster 



lOB, zero extinction would imply an age of -^20 Myr, an Ay 
of rjO.3 [E{B-V) « 0.10] intercepts the track at ~7-15 Myr, 
whereas an Ay ~ 0.8 [E{B — V) ~ 0.25] implies an age of 
~5 Myr. The extinction required for the combined point to 
cross the track at ~5-7 Myr is Ay m 1.6 [E{B -V) ^ 0.54]. 

Adopting a Galactic foreground reddening of Ay = 
1.64 \E(B - V) = 0.53] (following lOriglia et all 1200 ll : 
iRelano et al ]|2006l ) , we can immediately rule out ages requir- 
ing reddening values less than this. We therefore conclude 
that lOA must have an age of 5-7 Myr (with a reddening in 
addition to the Galactic foreground level of Ay ~ 0.0-0.4), 
and lOB an age of ^ 5 Myr (consistent with zero redden- 
ing after foreground correction). For the colour of the com- 
bined point to match these age determinations, the redden- 
ing value required is consistent with purely Galactic effects. 

A further test to obtain an estimate of the cluster age 
is to compare the equivalent width of the Ha emission line 
to the predictions of evolutionary synthesis models. Since 
we only have a spectrum of the combined light from both 
the clusters, we can only put an upper limit on the flux- 
weighted age. Measuring an equivalent width of Ho? = 91 ± 
3 A and using the same Sb99 model as used above {Z = 
0.004, Kroupa IMF), this gives an upper limit to the age of 
~6.5 Myr, and is consistent with the values derived above. 

Now we have an estimated age, we can derive an esti- 
mate for the mass by comparing the absolute magnitude. 
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3.1 Cluster Spectrum 
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Figure 3. Colour-colour diagram sliowing (solid line) Star- 
BURST99 instantaneous burst models for Z = 0.004. Crosses mark 
age points on the evolutionary curve from 5-7 Myr; some signif- 
icant age steps are labelled (units of Myr). The arrow represents 
the extinction vector for Apgssw = 1-0. Solid squares show the 
colours of the two clusters (cluster lOA and lOB); the solid circle 
shows the colour of the integrated flux of the two. 



.^Fsssw, to evolutionary synthesis models. As mentioned 
above, we cannot measure an absolute magnitude of the in- 
dividual sub-clusters, but using an aperture large enough 
to include both, we find a combined Mpsssw.ciusterio = 
—9.61 mag. Converting to Johnson magnitudes, we find 
Afi/,ciusterio ~ —9.2 mag, and comparing our measurements 
to a Sb99 6 Myr, 1 x 10® M0 model results in a predicted 
mass of 2-5 X 10^ M0 (depending on the IMF formulation 
used) for the combined mass of both sub-clusters. Since there 
are uncertainties involved in transforming from the HST 
to Johnson magnitude system, we compared our measured 
Afpss sw to the absolute magnit udes predicted by the equiv- 
alent iBruzual &: CharlotI (|2003l ) model (for which results in 
HST filters are given). The derived mass is 7-13 X 10^ Mq 
(depending on IMF), which is in general agreement with the 
Sb99 result. We therefore adopt the value of 7 ± 5 x 10^ M© 
for the combined mass of the two subclusters, meaning that 
neither of them can be defined as a super star cluster in the 
gen erally accepted definiti on . 

iBuckalew et al. I (|2000l ) conducted a study of NGC 1569 
using HST WFPC2 narrow-band imaging of the He 11 A4686 
line. They discovered fifteen Hell emitting sources across 
the galaxy, of which five were associated with stellar clus- 
ters including cluster 10 (specifically only the south-western 
component; lOB). They attribute the emis sion as nebular in 
origm (after iKobulnickv fc Skillmanlll997l . who found wide- 
spread nebular Hell emission in this region). However Hell 
is also strongly emitted in the atmosphere s of Wolf-Rayet 
(WK) stars of both WC and WN type. iBuckalew et all 
l|2000l ) comment that if the emission they detect in cluster 
10 was attributed fully to WR stars, it would correspond 
to the eq uivalent of three WN L-type stars (using the cali- 
bration of IVacca fc Contilll992l ). WR stars exist in clusters 
between ages of approximately 3-5 Myr, so in light of our 
age estimates for the two cluster components derived above, 
a WR origin for the He 11 emission in cluster lOB is perfectly 
reasonable. 



To extract the spectrum of cluster 10, we summed 29 spaxels 
(covering an area with an equivalent diameter to the FWHM 
of the cluster measured from the IFU continuum image) 
centred on the continuum source as shown on the right-hand 
plot of Fig. [5] The resulting spectrum is shown in Fig. 3] 

The spectrum is dominated by strong emission lines 
arising from the nebula surrounding the cluster, but also 
contains a broad emission feature that we identify with 
CivAA5801,5812 arising from WC-type Wolf-Rayet stars. 
We also detect a large number of weak stellar absorption 
lines, diffuse interstellar bands (DIBs) and interstellar ab- 
sorption lines. The detection of a red WR-bump is not un- 
surprising considering the age we find for cluster lOB, and 
strengthens the argument for a WR ori gin for the He 11 emis- 
sion detected by IBuckalew et al] (|2000l ). Further evidence of 
the young nature of cluster 10 comes from the high ratio 
of [Olll] to H/3 (see Fig. |4l lower line), and the presence of 
strong [S III] emission, indicating a high level of ionization 
in the surrounding Hll region, which can only arise in the 
presence of young O stars. 

Simultaneously fitting the WR feature with the super- 
imposed absorption lines, we measure the bump to have a 
FWHM of ~130 A, a total luminosity of 2.5 x 10^^ erg s"^ 
(assuming a distance of 2.2 Mpc and Ay = 1.64), and 
an equivalent width of 12 ± 2 A. This measured FWHM 
is unusually broad compared to can onical measurements 
(|Crowther. De Marco fc Barlow! 1 19981 ) , and is likely to re- 
sult from inaccuracies in the fit due t o contamination from 
the su perposed DIB absorption lines. ICrowther fc Hadfieldl 
(|2006l ) present measurements of WR bump fiuxes for a num- 
ber of WC4 type stars in the LMC (the closest galaxy con- 
taining a significant population of WR stars at low metal- 
licity - Z = 0.008) which can be used to determine the 
equivalent number of stars in our observation. They find 
one LMC WC4 star to have a red WR-bump luminosity of 
3.3±1.6 x 10"^® erg s~^, thus the luminosity equivalent num- 
ber of WC4 stars in our observed bump is estimated to be 
10 ± 5, once all the un certainties have been taken int o ac- 
count. For comparison. IConzalez-Delgado et al.l (| 19971 ) find 
a luminosity equivalent of 20 ± 10 WNL stars in SSC A by 
measuring the flux in the blue WR-bump. 

By fitting Gaussian profiles to all the other emission 
and absorption lines present in the spectrum, we mea- 
sured the central velocity and FWHM of each, and present 
them in Table [2] The average radial velocity of the nebular 
emission lines is —82 ± 4 kms~^. Since our cluster spec- 
trum only contains blended photospheric absorption lines, 
it is not possible to measure an accurate velocity for the 
stellar component. The evolutionary synthesis models of 
iBruzual fc CharlotI (|2003h provide a set of spectral energy 
distributions (SEDs) calculated for similar input ranges as 
the Sb99 models, and although being less accurate at pre- 
dicting very young ages, their SEDs are derived from actual 
observed spectral libraries rather than model stellar atmo- 
spheres. Their inclusion of observed stellar spectra at high 
resolution therefore make them ideal for comparing with 
our spectrum to derive a cluster velocity. Using a 0.25 Zq, 
5 Myr instantaneous burst model in the wavelength range 
5100-5600 A (where many faint stellar absorption lines are 
present), we find a match at —85 ± 5 kms~^. A good cor- 
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Figure 4. Plot of the extracted spectrum of cluster 10 (smoothed by 5 A). All identified absorption and emission features are labelled. 
The lower line corresponds to a x80 reduction of the main plot to illustrate the the spectrum with a full flux range. Notice the strength 
of the [Oin]A5007 emission compared to the hydrogen recombination lines. 



relation at this velocity is also found at the position of the 
Ca l+Fe blend at 6494 A, which is th ought to be photo- 
spheric in origin (|Heckman et al.lll995l ). As well as stellar 
absorption lines, we detect the interstellar (IS) abs orption 
lines of Nal AA5890,5896 and a number of DIBs (|Herbiel 
1 19951 : [ Heckman fc Lehnerd [20001 ). Their velocities are given 
in Table (2] and are consistent with being Galactic in origin 
(average velocity — 40±5 kms^^), although the Nal profiles 
also show blue wings presumably associated with IS absorp- 
tion within NGC 1569. 

Determining the proper heliocentric systemic velocity, 
Days, of NGC 1569 has proven to be a difficult task due to 
its irregular and disturbed morph ology. The fir st study of 
NGC 1569's H l gas distribution bv lReakesI (|l980|}_ concluded 
a Usva = —77 kms~^. CO observ ations by Travlor et al.l 



199£) foun d the same valu e, whils t Fibmit a. Ohta . 



19941 ) and iHeckman erahl 1^199^ ) found an Ha 



"Saitd 

^ of 

—90 kms~^ and —68 kms~^ respectively from their long- 
slit data. A more recent, high r esolution H I s tudy of the 
NGC 1569 system was made by iMiihle et all (120051 ). who 
found the peak of the Hi profile at the position of cluster 
10 to have a velocity of —80 kms~^ (S. Miihle, private com- 
munication). Thus, since our measurements of the ionized 



gas and cluster velocities are also consistent with this value, 
we adopt a systemic velocity for the cluster 10 region of 
-80 kms"\ 



fsys.lO 



A check on the value of reddening for the combined light 
of both cluster components can be made by comparing the 
obser ved Hq/H/3 fiux ra t io wit h the intrinsic case B values 
from iHummer fc StorevI |l983). We derive an absorption- 
corrected reddening of E{B — V) — 0.54 ± 0.15 mag from 
Gaussian fits to the Ha emission and H/3 emission and ab- 
sorption lines. This value is consistent with both the adopted 
value of Galactic foreground reddening and with the ex- 
tinctions derived from the colour-colour plot analysis given 
above. We therefore conclude that cluster 10 only suffers 
from Galactic foreground reddening, i.e. very little of its 
light is attenuated by material in NGC 1569 itself. The flux 
ratio of single Gaussian fits to the [S ll]AA6717,6731 emis- 
sion lines gives a value of 1.25lJ!;5;i hence an electron den- 
sity, Ua — 175^20 cm~^ (assuming an electron temperature, 
Tc = 10" K). These resuhs are summarised in Tabled 
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Table 2. Absorption and emission line measurements for clus- 
ter 10 (all lines are emission unless otherwise stated). Widths arc 
corrected for instrumental contribution and velocities are quoted 
in the heliocentric frame of reference. Typical errors on the veloc- 
ities of the emission lines are ~5— 10 kms~^, and on the widths 
~0.5-5 kms~^. Errors for the absorption lines are listed individ- 
ually. 



Table 4. Wavelength limits used in fitting each line 



Spectral Line Wavelength (A) 



Line 


-^air 


Velocity 


FWHM 




(A) 


(kms-l) 


(kms-i) 




4Sfi1 S'^ 


—84.4 


50.7 




^OU -L . OO 


—55 it 40 




T-To T aVid ^ 


4921.93 


—84 ± 8 




[Oiii] 


4958.92 


-79.0 


44.7 


[O III] 


5006.84 


-75.1 


31.6 


Hei'= 


5015.68 


-74.7 


42.3 


Mgl abs 


5167.32 


-77 ±7 




Mgl abs 


5172.68 


-89 ±6 




DIB 


5780.5 


-40.9 ±10 


160 ± 35 


DIB 


5797.0 


-47.5 ±20 


70 ± 25 


He I 


5875.67 


-82.7 


33.9 


Nal abs 


5889.95 


-38 ±8 




Nal abs 


5895.92 


-34 ±8 




DIB 


6283.9 


-41.0 ±10 


110 ± 20 


[Oi] 


6300.30 


-84.9 


32.8 


[Siii] 


6312.10 


-82.0 


50.5 


[Oi]" 


6363.78 


-80.1 




[Nil] 


6548.03 


-81.4 


28.3 


Ha* 


6562.82 


-89.1 


38.8 


[Nil] 


6583.41 


-82.8 


44.5 


He I 


6678.15 


-81.6 


44.1 


[Sii] 


6716.47 


-85.1 


40.8 


[Sii] 


6730.85 


-84.4 


46.5 



" Line too weak or confused to measure an 
accurate FWHM 
' Narrow component (CI) only 
Blended with [O ill] 



Table 3. Observed and derived properties for cluster 10. 



Parameter 



Cluster component 
lOA lOB 



Adopted systemic velocity 


-80 


(km s~^) 




Distance (Mpc) 


2.2 ±0.6 + 


Coordinates (J2000) 


04^ 30™ 47? 16 04^ 30™ 47! 12 
+64° 51' 00'.'8 +64° 51' 00'.'6 


(F555W-F814W) 


1.12 0.42 


(F330W-F439W) 


-0.27 -0.90 


Half-light radius, R^ff (pc) 


0.88 + 0.05 0.60 ±0.15 


Total E{B - V) 


0.55 ±0.15 


Electron density, no (cm~^) 


1 7^ + 18 


Age (Myr) 


5-7 ^5 


A^F555W (mag) 


-9.61 


Mass (Mq) 


7 ± 5 X 10^ t 


No. of equiv. WC4 stars 


10 ±5 


Radial velocity (kms~^) 


-85 ± 5 



t [Israeli l|l988h 

t Depending on whether a Salpeter or Kroupa IMF is used 



H/3 


4840 


-4880 


[O iii]A5007 


4940 


-5020 


[Oi]A6300 


6250 


-6308 


[N ii]A6583 


6570 


-6595 


Ha 


6550 


-6575 


[S ii]A6717 


6700 


-6724 


[Sn]A6731 


6720 


-6750 



4 PROPERTIES OF THE IONIZED GAS 

We now discuss the IFU data relating to the ionized gas in 
the environment of cluster 10. The signal-to-noise and spec- 
tral resolution of these data are sufhciently high to resolve 
multiple Gaussian components to each emission line in each 
of the 500 spectra across the field-of-view. Before discussing 
our findings, we first describe the methods used to fit the 
spectral line profiles and visualise the results. 



4.1 Decomposing the line profiles 

The automated fitting of multiple profiles to 500 spectra 
each containing eight emission lines has proved challeng- 
ing. For this task we have adopted an IDL-based, general- 
purpose cur ve-fitting package, called PAN (Peak ANalysis; 
iDimeollioOsI ). It was written with the aim of allowing the 
user to visually interact with the program through the use 
of IDL's sophisticated 'widget-based' toolkits, and uses the 
Levenberg-Marquardt technique to solve the least-squares 
minimisation problem. The primary features that caused us 
to favour this program over any other was its ability to 
read in multiple spectra at once in an array format, and 
its method of determining the initial guess parameters of 
a model profile (namely allowing the user to visually and 
interactively specify amplitude, position and width). Orig- 
inally written to analyse data from neutron-scattering ex- 
periments, the program has required some modification to 
suit our needs. We have had to add modules to convert our 
FITS format data into a format that PAN can read; we have 
also had to modify the way it outputs the fit results and 
values to something that can easily be read by a plotting 
package. 

To obtain good, reliable fits, the spectra read into PAN 
had to contain accurate associated error arrays. These error 
arrays were calculated by taking the square-root of the raw 
counts in the pre flux-calibrated and throughput-corrected 
data, converting this to a percentage of the total raw counts, 
and applying this to the final reduced and calibrated data in 
order to obtain the uncertainty equivalent to the Poissonian 
noise. This step was included when converting from FITS to 
PAN input format. By comparing our error estimates on fits 
to telluric features (see below) to simple Poissonian statis- 
tics, we find this to be a realistic measure of the cumulative 
uncertainties. 

Each line in each of the 500 spectra was fitted us- 
ing a single, double and triple Gaussian component ini- 
tial guess, where line fluxes were constrained to be positive 
and widths to be greater than the instrumental contribu- 
tion (measured from wavelength calibrated arc exposures 
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to be 74±5 kms~^ at [Olll] and 59±2 kms"^at Ha; Sec- 
tion [2]3]). The wavelength limits used in fitting each line are 
given in Table U) In general, the shape of the emission line 
profiles is a convolution of a bright, narrow component over- 
lying a faint, broad component (most obvi ous in high S/N 
lines, e.g. Ha). This is consistent with what iHeckman et al.l 
j 19951 ) found for the central regions. For the double compo- 
nent fits, the initial guess was always made with the first 
Gaussian as the narrow component (hereafter referred to as 
CI), and the second Gaussian as the broader component 
(hereafter C2). For the triple component fit, the additional 
component (C3) was specified with a initial guess assigning 
it to a supplemental narrow line at the same wavelength as 
the main narrow line. This consistent approach helped limit 
the confusion that might arise during analysis as to which 
Gaussian fit belonged to which component of the line, as 
well as aiding the x'^ minimisation process. 

We used the statistical F-test to determine how many 
Gaussian components best fit an observed profile. The F- 
distribution is a ratio of two distributions, denoted by 
the degrees of freedom for the numerator and the de- 
nominator X . For statistically comparing the quality of two 
fits, this function allows one to calculate the significance of 
a variance (x^) increase that is associated with a given con- 
fidence level, for a given number of degrees of freedom. The 
test will output the minimal increase of the x^ ratio that 
would be required at the given confidence limit for deciding 
that the two fits are different. If the x^ ratio is higher than 
this critical value, the fits are considered statistically distin- 
guishable. The output files for the single, double and triple 
Gaussian models were passed through a script which applied 
the F-test to extract the results for the most suitable case 
for each spectral line. Fig. [S]shows two example Ha line pro- 
files, each from a different spaxel, with their corresponding 
Gaussian fits and residual plots shown. PAN calculates the 
residuals, ri, using the following formula 



where ai are the uncertainties on yf^'*. The upper exam- 
ple represents an instance where two Gaussian components 
are needed to fit the profile — with the inclusion of a second 
Gaussian the fit x^ is dramatically reduced and the resid- 
uals fall within the ±1 guidelines. The lower plots show an 
example of a profile which is best fit with a triple Gaussian 
model. The addition of a second narrow, bright component 
as well as the broad, faint component causes a significant 
reduction in the x^ and again the residuals fall within the 
±1 guidelines. 

We also applied a number of physical tests to further 
improve the accuracy of the results. Firstly, for a fit to be 
accepted, the measured FWHM had to be greater than the 
associated error on the FWHM result (a common symptom 
of a bad fit), and secondly the FWHM had to be smaller 
than 15 A (to guard against fitting continuum or spurious 
features). We also applied some additional filtering in order 
to ensure the first Gaussian component always corresponds 
to our definition of the first component (i.e. if fiux(Cl) < 
flux(C3) AND fwhm(Cl) < fwhm(C3) then swap 03 with 
CI). This only applies when the line contains three com- 
ponents and is sufficiently well resolved to reliably fit them 
aU. 



The bright [Olll]A5007 line also exhibits the general 
line profile shape described above (narrow bright compo- 
nent with underlying fainter broad component), but unfor- 
tunately due to an instrumental artefact present at low lev- 
els in the blue wing of the line, we have not been able to 
characterise the shape of the secondary components accu- 
rately. The flux, width and velocity maps of the main bright 
[O ni]A5007 component (CI) look very similar to the equiva- 
lent Ha maps, so we do not show them here. Measurements 
of the [O lll]A4959 hue also result in similar flux, FWHM and 
velocity distributions, but since it is a weaker line we cannot 
fit the low intensity broad component with any confidence. 



4.2 How well can we fit each component? 

The x^ minimisation algorithm used by PAN is very sensitive 
to the uncertainties associated with the input spectra. Be- 
cause of this, we paid special attention to creating accurate 
Poissonian error arrays from the raw un-calibrated spectra, 
that we then associated with the fully reduced spectra for 
the purposes of fitting (see Section|4]T}. We made a number 
of tests comparing the results of Gaussian fits to high S/N 
arc lines made with PAN with those output by other line 
fitting routines (such as elf within the STARLINK dipso 
program). We found that despite differences in the minimi- 
sation techniques employed, results and uncertainties quoted 
were very similar, and that PAN produced more consistent 
values across the whole wavelength range. We are therefore 
confident with the quality of results derived by PAN. 

For a consistency check on the variability of the 
instrumental profile, we performed Gaussian fits to the 
unresolved [Ol]A5577 telluric emission line on pre-sky- 
subtracted spectra. The spatial distribution of the fit 
components clearly show no systematic variability across 
the IFU field indicating that the line profile shape is 
very stable. We also find that the telluric line is well 
fit by a single Gaussian, and shows no evidence for 
broad wings. The [O l] line fluxes have a standard devia- 
tion, s.d. = 1.5 X 10^^ erg s"^ cm"^ arcsec"^, the FWHM 
s.d. — 4.5 kms^^ and radial velocity s.d. = 1.9 kms^^. 

We now turn to the uncertainty estimates on each 
line component's Gaussian properties (line centre, flux and 
width). Unfortunately, the errors that PAN quotes on its fit 
results are derived simply from the formal errors on the 
minimisation, and we have found these to be an under- 
estimate of the true uncertainties. More realistic errors have 
been estimated through visual inspection of the fit quality 
taking into account noise in the continuum (i.e. S/N of line) 
and the number of fitted components. For the Ha line flux, 
the percentage error varies between 0.5-10 per cent for CI 
(for high-low S/N lines), 8-15 per cent for C2 and 10-80 
per cent for C3. The addition of multiple Gaussian compo- 
nents always increases uncertainties; we estimate that where 
we can fit a third component, the errors in the flux of CI 
and C2 increase by 5-10 per cent. For estimation of the 
FWHM and radial velocity errors, we have compared fits to 
the two lines of the [O lll]AA4959,5007 doublet where possi- 
ble. We find the FWHM errors vary between 0.25-2 kms~^ 
for high-low S/N CI lines respectively, 4-15 kms"'^ for C2 
and 15-20 kms"'^ for C3. For the radial velocity, the error in 
CI varies between 0.1-3 kms~^, 2-5 kms~^ for C2 and 10- 
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Figure 5. Examples of Ha line fits made by PAN comparing a single and double Gaussian fit for spaxcl 114 and a double and triple fit 
for spaxel 449 (j/-axis in arbitrary fiux units, but note the difference in scale). The black solid line is the observed data and the dashed 
lines represent the individual Gaussian components. Below each graph is a plot of the residuals (in units of cr) with dashed guidelines 
plotted at ±1 (see text for an explanation of how the residuals are calculated). By comparing the ratios of the values, the significance 
of a higher number of fit components can be tested. 



30 kms~^ for C3. We estimate that the addition of a third 
component (where required) increases the errors in the CI 
and C2 FWHM by ~5 kms"^ and the central velocity by 
^3 kms^^. 



4.3 Visualisation of the Data Products 

A major hurdle in the analysis of IFS is the visualisation 
of the data products. To this end, we have developed a 
PERL/PGPLOT based plotting and analysis package nick- 
named 'daisy' that has enabled us to map out the results 



that we have obtained from our dataset. We have already 
shown one such visualisation in Fig. [5] The program allows 
the user to correct for DAR (currently not implemented), 
correct the measured FWHMs for the intrinsic instrumen- 
tal profile (via simple FWHM^ subraction), deredden the 
data using the intrinsic Ha/H/? ratio method, convert wave- 
lengths and FWHM into velocities, create a number of ratio 
combinations, and then to visualise the results. There are 
options to allow the user to select the colour range (linear 
or log) for each plot, and to overplot using contours. Once 
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fully developed, our aim is to release this package to the 
communitjrl. 



5 EMISSION LINE MAPS 

Figs. E] to [8] show maps of the fitted Gaussian properties 
(flux, FWHM and radial velocity) for each ffa line compo- 
nent (we show maps for measurements of the Hq line only 
for the reasons given above in Section I4.1|l over the whole 
field. The grids represent that of our interpolated, cubed 
data, not the original GMOS configuration, and the posi- 
tion of the cluster is overlaid with flux contours on selected 
maps to aid description and interpretation. As described in 
Section [4.11 in general we find a bright, narrow component 
overlaid on fainter, broad emission. This Ha broad compo- 
nent reaches a corrected FWIfM of ~300 kms~^ in places. 
In a number of spaxels we have been able to fit a third com- 
ponent (C3), which has approximately the same width as CI 
implying that it may be associated with expanding shell- or 
bubble-like structures in the ionized gas. 



5.1 Ha Maps 

5.1.1 Component 1 (CI) 

From Fig. [S] it is immediately obvious that CI and C2 do 
not originate in the same place, i.e. the flux maps show strik- 
ingly different spatial distributions. The two main knots of 
emission in CI, located near the middle and south-west of 
the field (south of the cluster), correspond to the pattern of 
diffuse emission seen in both the WFPC2 F656N (Fig. [T]) 
and ACS F555W (Fig.O images. These knots lie just to the 
south of cluster 10 and mark the easternmost exten t of the 
large ionized complex seen in Fig. [T] (H ll region 2; IWalled 
1 19911 ). The width of CI shown in Fig. |7] (left) ranges from a 
FWHM of ~30-40 kms~^ (corrected for instrumental broad- 
ening) towards the eastern side of the field, to ~60 kms^^ 
seen in a region arcing round to the west of cluster 10. This 
is coincident with the strip in-between two linear filaments 
seen to the north-west of cluster 10 in Fig. [T] and may be 
indicative of an influence from the two SNRs located to 
the north of the FoV. A gradient in FWHM from ~35 to 
45 kms~^ is centred on the position of the bright Ha knot 
in the centre of the field, with broader widths seen on the 
side facing cluster 10, but the line width stays constant to 
within ~5 kms~^ over the extent of the cluster. As shown 
by the CI radial velocity map (Fig. [5] left), we see a re- 
gion of material moving at ~— 10 kms~^ in the centre (just 
to the south of the cluster) surrounded by a ring of mate- 
rial with a velocity decreasing down to the systemic value. 
There is also an area of gas just to the north moving with 
a radial speed of ~-|-7 kms~^. It is worth noting that these 
regions of negative and positive velocities are located either 
side of the position of the cluster at approximately equal 
distances. The fastest radially approaching CI gas is seen 
in the south-west, moving at 10-15 kms~^ (coincident with 
the small dense knot seen in the HST Ha image), and the 

^ interested users should contact Katrina Exter, kexter@stsci.edu 



fastest receding gas is found along the eastern edge, increas- 
ing to ~15-20 kms~^ (off the colour scale for this plot) in 
the far north-east. It should be borne in mind that the total 
CI radial velocity range observed across the whole field is 
only ~35 kms^^. 

5.1.2 Component 2 (C2) 

The two small bright knots seen in C2 (Fig. |B] centre) to- 
wards the north of the field (north-east of cluster 10) do 
not correspond to anything seen in the HST images. This 
is understandable since the most intense emission in CI is 
an order of magnitude brighter than that in C2. The dis- 
tribution of C2 FWHMs (Fig. [T] centre) is almost precisely 
anti-correlated with the line flux (Fig. [§1 centre) , where the 
brightest features exhibit the narrowest line width (115- 
125 kms~^). The left-hand panel of Fig. [9] better illustrates 
this by showing contours of C2 flux overlaid on the C2 
FWHM map. The broadest C2 lines (~300 kms~^) are seen 
in the south-centre extending to the south-west of the fleld, 
corresponding well with the position of the brightest knots 
seen in CI (Fig. [6] left) and the WFPC2 image. This again 
is illustrated clearly in Fig. [9] (right-hand panel), where con- 
tours of CI flux are overlaid on the C2 FWHM map. The C2 
radial velocity field (Fig. [S] centre) covers the same overall 
range as the CI gas; material approaching us at velocities of 
~15 kms~^ (shown in blue) is located near the northern face 
of the bright CI knot in the centre of the field, in-between 
the negative and positive velocity regions in CI, and almost 
coincident with cluster 10. This area of approaching gas also 
extends towards the west, reaching a maximum approach- 
ing velocity of 20 kms~^. The fastest receding C2 velocities 
(red) are found along the eastern edge of the field at ~20- 
25 kms~^ relative to the average. The spatial distribution 
of this region of redshifted gas mimics that of the redshifted 
CI gas, but extends further to the west. Note the distribu- 
tion of C2 radial velocities does not correlate with either of 
the CI or C2 flux maps, however we do tentatively point out 
a possible radial velocity difference of ~25 kms~^ between 
the two bright C2 knots seen in the flux map (Fig. |6l cen- 
tre), and a spatial coincidence with the position of VLA-16 
(a possible extended low surface brightness SNR). 

5.1.3 Component 3 (C3) 

A third component to the Ha line profile fit was required 
in a small region located in the north-east of the field. The 
flux map (Fig. [B] right) shows it is brightest in the far north- 
eastern corner and has a FWHM range (Fig.[7l right) similar 
to CI. The radial velocity map (Fig. [SI right) indicates that 
the C3 gas has a larger velocity towards our line-of-sight 
than any of the CI or C2 gas (up to ~70 kms~^), but the 
associated uncertainties are large (see Section 14. 2p . Since 
our only detection of this third component is very near the 
edge of the fleld, it is difficult to estimate its extent or role 
from just these maps. Perhaps more interesting is the non- 
detection of a third component in locations where signatures 
of expanding structures might be expected from inspection 
of Fig. [T] We discuss the significance of this component in 
context with the whole central region and wind-blown cavity 
in Paper II. 
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Figure 6. Left: Flux map in Ha component 1 (CI; bright, narrow); centre: flux map in Ha component 2 (C2; faint, broad); right: 
flux map in Ho component 3 (C3; narrow, second peak). Non-detections are represented as hatched spaxels, the x and y scales are in 
arcseconds offset from the lower-left spaxel, and a scale bar is given for each plot in units of 10~^^ erg s~^ cm~'^ arcsec"^. North is up 
and east is left. The higher levels in the continuum image are plotted as contours on the left-hand plot, identifying the position of cluster 
10. 
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Figure 7. Left: FWHM map in Ha CI; centre: FWHM map in Ha C2; right: FWHM map in Ho C3. A scale bar is given for each plot 
in kms^^ units corrected for instrumental broadening. North is up and east is left, and the position of cluster 10 is identifled on the 
left-hand plot. 



5.2 Flux ratio maps 

5.2.1 Electron density 

Fig. [To] shows a plot of the spatial variation of electron den- 
sity, rio, as derived from the ratio between single Gaussian 
fits to the [S ll]AA6713, 6731 lines, assuming an electron tem- 
perature, Tc, of 10* K (the S/N of these lines is not high 
enough to fit multiple line components). There are two dis- 
tinct regions of higher density seen in the map: in the south- 
west of the field, Uc reaches ~250 cm~^ and is spatially co- 
incident with the south-western bright knot seen in Ha CI 
(Fig. [6] left), and the small distinct knot seen in the HST 



F656N image. A second high density region is coincident 
with the north-eastern part of the central bright Ha CI knot 
and the position of cluster 10, and peaks at ~270 cm^^. We 
derived the electron density of the Hll region surrounding 
cluster 10 in Section [3] and found = 175^20 cm^"^. Our 
peak value of 270 cm~^ is entirely consistent with this value 
since the spectrum for cluster 10 was created by summing 
29 spaxels centred on this central density peak, and hence 
includes some dilution from the surroundin g lower density 
gas. T hese values are also consistent with iHeckman et al.l 
(| 19951 ). who find densities of up to 200 ± 100 cm"^ in the 
central 15" of NGC 1569. Only a small proportion of the gas 
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Figure 8. Left: radial velocity map in Ha CI; centre: radial velocity map in Ho C2; right: radial velocity map in Ha C3. A scale bar is 
given for each plot in units of kms~^ (heliocentric) relative to the systemic velocity of the region (—80 kms~^), where zero is marked 
with a line. The position of cluster 10 is identified on the left and central plots for comparison. 




Figure 9. Lefl.: contours of C2 flux overlaid on the C2 FWHM map. Solid contours represent levels above the mean of the range, and 
dashed below the mean. Notice the precise correlation between areas of high flux and low FWHM in the north. Right: contours of CI 
flux overlaid on the C2 FWHM map. Here the correlation between the bright Ho knots seen in CI and the broadest C2 lines can be seen 
clearly. 



is at or below the low density limit (~100 cm ^: lOsterbrockl 
1 19891 . shown by the single contour level on the map). 



5.2.2 Nebular diagnostics 

Our wavelength range gives us access to a number of im- 
portant optical diagnostic lines that can be used to con- 
strain the properties of the emitting gas. The flux ratio of 
[Olll]/H/3 is a good reddening- free indicator of the mean 
level of ionization (radiation fleld strength) and tempera- 
ture of the gas, whereas [Sll]/Ha or [Nll]/Ha are indica- 



tors of the num ber of ionizations per unit v o lume (ioniza- 
tion parameter ; IVeifleux fc Osterbrockl [l987l : iDopita et al.l 
A greyscale representation of the observed CI flux 
ratio of [O lll]A5007/H/3 overlaid with contours of the 
[Sll](A6717-hA6731)/Ha CI flux ratio is shown in Fig. [TT] 
The ratio of log([S ll]/Ha) is lowest ( — 1.7; shown by dashed 
contours) where the flux in Ha CI is highest (see Fig. |6l 
left), i.e. where the central bright knot is located. The peak 
of the [Olll]/H/3 ratio is also near to the central bright 
Ha knot, but offset by ~4'.'5 to the south-east. Numeri- 
cally, the field average ratio of log([S ii]/Ha) = —1.29 ± 
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Figure 10. Map of the electron density, no, calculated from 
the ratio of single Gaussian fits to the the [Sll]A6716, 6731 lines. 
The scale is in units of cm~^. A single contour level marks the 
100 cm~^ low density limit. 

0.17, log([Nn]/Ha) = -1.62 ± 0.16, and log([0 iii]/H/3) 
= 0.87 ± 0.05 (one a errors), showing that the variations 
in the ratios are really only small. Plotting these values 
on the st andard diagnostic diag r ams o f the type first pro- 
posed by IVeilleux fc Osterbrockl (|l987l ) , shows that all the 
points lie tightly clustered in the upper-left of the starburst- 
Hll region distribution. Comparing t hese values with the 
Mappings-Ill photoionization models of lDopita et al.l (|2000l ) 
for an instantane ous burst of star-form ation whose SED 
is given by Sb99 (|Leitherer et al.lll999l ), the gas is of low 
metallicity (~0.2 Zq) and has a high ionization param eter 
[q ~ 10®), as expected. iBuckalew fc Kobulnickvl (|2006h use 
//Sr/WFPC2 narrow-band imaging to calculate the ratios 
of [Olll]/H/3 and [Sii]/Hq:, a nd use the theoretica l 'maxi- 
mum starburst line' derived bv lKewlev erall l|200lh to set a 
threshold to determine the difference between photoionized 
and non-photoionized emission. Using this criterion on our 
de-reddened line ratios, we find a total of 118 spaxels con- 
taining possible non-photoionized emission, the location of 
which are plotted on Fig. llllwith white cros ses. We note here 
that the threshold derived bv lKewlev et al.l is a maximum in 
the sense that no pure starburst, regardless of metallicity (in 
the range 0.1-3 Zq) or ionization parameter, could produce 
emission above this line. A lower metallicity, such as found 
in NGC 1569, decreases the level of this theoretical limit, but 
since the model predictions are quite uncertain, understand- 
ing how many more points may have a non-photoionized 
contribution to their emission is non-trivial. For this reason, 
we have employed the more robust 'maximum' threshold. 

The distrib ution shown in the map corres ponds well 
to the maps of iBuckalew fc Kobulnickvl (|2006l ). with the 
majority of points lying in the south-east of the field. 
Since SSC A is located only 7" away in this direction, 
IBuckalew fc Kobulnick"^ (|2006l l associate these points with 
an 'arc' of non-photoionized gas resulting from a possible 
wind-wind interaction between the two clusters. 



arcsec 

Figure 11. Flux ratio of log([0 in] A5007/H/3) CI 
(grcyscale; scale bar shown) with contours of the ratio 
log([Sll](A6717-(-A6731)/Ha) CI (12 linearly spaced levels 
from —1.7 to —0.9). Spaxels marked with a cross indicate 
non-p hotoionized emission according to the maximum starburst 
line of lKewlev et al.l l|20o3). 

5.2.3 Reddening map 

De-reddening with our v isualisation and analysis package 
uses the Galactic law of iHowarthI (Il983|. Ry = 3.1) and 
the in trinsic Hi ratio of Ha/H/3 = 2.85 ( Storev fc Hummeij 
1 19951 . for gas at Te = 10* K and = 100 cm"^). A redden- 
ing map, computed using only CI of Ho? and H/3, is shown 
in Fig. 1121 The mean and standard deviation for c(H/3) for 
CI over the field-of-view is 1.01 ± 0.16 (compared to 1.15 
when summing the flux in all components for the H I lines) . 
In this, we are assuming most of the reddening c omes from a 
Galactic foregroun d screen (a valid assumption; Ilsrael|[l988l : 
lOriglia et al.ll200ll . Section The reddening distribution is 
consistent with the location of the central and south-western 
bright knots seen in CI, where these knots exhibit the high- 
est extinction [c(H/3) = 1.3 = E{B-V) = 0.8]. Interestingly, 
high values are also found in the north-west of the field, not 
correlating with any feature on any of the other maps. The 
lowest reddening [c(H/3) = 0.75 = E(B - V) = 0.5; equiv- 
alent to the Galactic foreground level; iRelano et al.l [20061 ] 
is found in the north-east of the field, coincident with the 
location of the group of spaxels containing C3. 



6 INTERPRETATIONS AND DISCUSSION 

At the distance of NGC 1569, the CMOS IFU field-of-view 
(FoV) covers ~50 x 35 pc with an 8 pc resolution, and is 
therefore ideal for studying the detailed relationship between 
young star clusters and their environment. To this end, we 
have observed a region near the centre of the galaxy covering 
one such young, bright star cluster and the surrounding ion- 
ized gas. Taking advantage of the spatially resolved nature of 
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Figure 12. A c{H/3) map calculated from the CI flux ratio of 
H/3/Ha. 

IFU observations, we extracted a spectrum of the cluster to 
derive its fundamental properties, and created well-defined, 
high S/N maps of the Gaussian properties of each observed 
Ha line component across the FoV. In the following sections 
we discuss our findings relating to the cluster itself and to 
the surrounding ionized gas. We also explore possible the- 
ories that might give rise to the observed broad emission 
profiles and discuss their consequences. 



6.1 Cluster 10 

As clearly seen in the high-resolution HST / ACS images 
(Fig. [2}, cluster 10 is composed of two individual clusters 
which we denote lOA and lOB. We measure their projected 
separation to be ~3.7 pc and their sizes to be ReS,a = 0.88 
and ReS,B = 0.60 pc, and a combined photometric mass of 
7±5 X 10'^ Mq. By comparing the relative photometry of 
each cluster we find them to have very similar ages: lOB is 
the youngest (<5 Myr) and lOA a little older (5-7 Myr). 
The combined spectrum of the two clusters extracted from 
our IFU data shows that one (or both) of these clusters has 
a significant population of WR stars. 

The unusually compact nature of these two clusters is 
reminiscent of R136 in the LMC. This equally young clus- 

ter has an R eH « 1.7 pc and mass of a few x 10'' Mq 

jHunter et al.l 1 19951 : iBosch et all bOOll ). R136 contains a 
large number of massive stars that power and maintain the 
~100 pc-diameter giant Hll region 30 Dor. 

Our age estimates of the two clusters show that their 
formation was approximately coeval with SSC A, and that 
vigourous star-formation has certainly occurred recently in 
NGC 1569. Cluster 10 is located within the bright H ii region 
No. 2 (|Walleijll99lt l which is found to e mit a powerful ther- 
mal radio signature (|Greve et al.ll2002l) and bright mid-IR 
forbidden-line and continuum emission (jTokura et al.ll2006h . 



These observations imply that this whole area is a large em- 
bedded star-formation compl ex being fuelled by the neigh- 
bouring molecular CO cloud (|Tavlor et al.|[l999l l. 

We do not see evidence for a compact individual Hll 
region surrounding the cluster in the Ha flux maps or the 
extinction map. This is probably because the winds and su- 
pernova events have blown away the surrounding gas. The 
electron density map, however, does show an enhancement 
at the location of cluster 10, indicating that perhaps there 
is some associated gas remaining. Our IFU FoV contains 
a great deal of ionized material that is both influenced by 
cluster 10 and by the relatively distant, but more massive, 
SSCs A and B. It is not clear if the source of ionization for 
Hll region 2 is solely cluster 10 or whether other starburst 
regions contribute. 



6.2 Ionized gas 

6.2.1 The turbulent ISM and the narrow component (CI) 

The integrated Hi velocity dispersion for the disc of NGC 
1569 is o- f» 15 kms"^ (FWHM 3 5 kms~M indicatiri g that 
the neutral ISM is very disturbed (jMiihle et al.ll2005h . pre- 
sumably due to the combined effects of stellar winds and su- 
pernovae. This value is consistent with the broad line widths 
also fou nd in Hi clouds associated wi th nearby giant Hll 
regions (|Viallefond. Allen fc Go"s3ll98ll ). We would thus ex- 
pect the Hll gas in NGC 1569 to have similar or higher 
velocity dispersions to the H I. 

In Fig.[l3]we plot the heliocentric radial velocity against 
line width for all identifled line components across the FoV, 
whilst in Fig. [14] we plot the line flux against line width. 
This latte r plot was used su ccessfull y as a dia gnostic tool by 
Mufioz- Tunon et al ] (| 19961 . see also iMartinez- Dclaa do et al.l 
2007i ) to examine the physical origin of the broadening mech- 
anisms in giant H ll regions. Plotting the line characteristics 
in this way highlights a number of trends within and between 
the components. 

Fig. [T3l shows that components CI and C3 (FWHM = 
35-65 kms~^) are easily distinguished from component C2 
(FWHM = 120-300 kms"^). Furthermore, in the bottom 
panel of Fig. 1141 a distinct ~35 kms~^ lower limit to the CI 
FWHM (which we henceforth refer to as FWHMo) can be 
seen. After correction for thermal broadening (~20 kms^^at 
T ~ 10* K), this becomes 29 kms"^ (or a ^ 12 kms"^), 
meaning the CI line velocity widths correspond to mildly 
supersonic speeds (for Tgas = 10* K, the sound-speed is 
~10 kms~'). We note that the CI line width does not cor- 
relate strongly with line flux, unlike the broad component 
(C2) which we discuss below. 

The red velocity tail of CI, seen clearly on Fig. 1131 as 
points lying to the right-hand side of the vertical dashed line 
at velocity = -1-10 kms~^, corresponds spatially to the ex- 
treme north-east of the field where we also see a blueshifted 
component (C3; see e.g. Fig. [5] lower panels). These C3 
points lie to the left-hand side of dashed line at velocity 
= — 20 kms~*. The velocities of these double-peaked emis- 
sion lines are distributed evenly to the red and blue of the 
cluster radial velocity (shown as a grey strip), providing the 
only signature of classical shell expansion in the field. Fur- 
thermore, the two components have similar fiux levels as 
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Figure 13. Ha line width (FWHM, corrected for instrumental 
broadening but not thermal broadening of ~20 kms~^) vs. cen- 
troid velocity (heliocentric, relative to the systemic velocity of the 
region, Vays = —80 kms~^). A representative error bar has been 
plotted for each of the three components on one selected data 
point. The shaded region represents the velocity (including the 
error) of cluster 10 (see Section Next to the right-hand axis 
are labels marking the three regimes in FWHM discussed in the 
text. CI (C3) points lying to the right (left) of the vertical dashed 
lines form the high- velocity tails referred to in the text. 
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Figure 14. Ha line width (FWHM, corrected for instrumental 
but not thermal broadening) vs. flux (arbitrary scale). The se- 
lected axis ranges show CI and C3 in the lower panel, and C2 
in the upper panel (note in particular the factor 10 difference 
in flux range). The distinct lower limit to the FWHM of CI, 
FWHMq Ri 35 kms~^, is marked. Encircled points are those 
belonging to the high-velocity tails referred to in the text and 
defined by the vertical dashed lines in Fig. 1131 



illustrated in Fig. [14] (see also Fig. [5] lower- right panel), 
where the corresponding high- velocity CI and C3 points are 
plotted with encircled symbols. 

Broad line widths have been observed in many 
young star-formin g regi o ns e.g. 30 Dor o r NGC 604 
\c\m fc KennicuttI Il994l : lYang et all Il996l ). A search 
through the literature uncovers a number of alternative the- 
ories that could account for the observed line- widths of CI, 
including electron scattering, multiple unresolved expanding 
shells integrated along the line-of-sight, and gravitational 
broadening through virial motions. 

Electron scattering can be immediately ruled out, 
as unphysical electron d ensities are requ ired to produce 
the observed widths (see IRov et al.|[l992l ). Multiple unre- 
solved expanding shells integrat ed along the line-of-sight 
was an explanation proposed by IChu fc Kennicutg (|l994h 
to explain their observed broad line-widths in the giant 
Hll reg ion 30 Dor. A follow up study of the s ame re- 
gion bv iMelnick. Tenorio-Tagle fc TerlevichI ljl999l ). but at 
very high spatial- and spectral-resolution (O'.'G seeing and 
0.067 Apix~^), allowed shells to be resolved down to their 
resolution limit of 0.1 3 pc. In their observatio ns, the broad 
line- widths found by IChu fc Kennicutg Il994l were resolved 
into a myriad discrete narrow lines representing individ- 
ual shell components, thus validating the conclusions of 
IChu fc KennicuttI (|l994h when observations are made at in- 
sufficient spatial / spectral-resolution. 

However an analyti cal model for this unresolved shells 
effe ct, first developed by DvsonI (Il979l) and analysed i n deta il 
by iTenorio-Tagle. Munoz-Tunon fc Cid-Fernanded l|l996t ). 
predicts that if the age distribution of the stars powering 
the expanding shells is greater than the timescale for gravi- 
tational collapse of the parent cloud, the resulting line pro- 
files are always flat-topped. This is something we do not 
see. Since the spatial and spectral resolution of our obser- 
vations is ~9 pc and 0.34 Apix~^, it is conceivable that a 
contribution to the narrow-line width is due to unresolved 
expanding shells. However, since the point at which coher- 
ent shells break up is when their expansion velocity equals 
the ambient ISM turbulent velocity (FWHMo), unresolved 
expanding shells cannot contribute significantly to the ob- 
served widths of CI; 

Furthermore, iTenorio-Tagle et all (|l996l ) find that to 
have a collection of unresolved shells capable of producing 
a supersonic line-width over a large range of densities, im- 
plies a great many low-energy winds with properties similar 
to those typical for low-mass stars. If multiple shells are re- 
sponsible for the observed broadening in our case, then the 
level of broadening would be expected to scale to a certain 
degree with gas column density (a large volume of gas could 
contain shells with a larger range of velocities), but this is 
not observed. The broadest CI lines are seen in the north- 
west of the field, away from the bright knots in the field 
centre. 

Gravitational broadening has been suggested by a num- 
ber of authors as an ex planation for th e observed line- widths 
in giant Hll regions (jM clnick' '1977"; 'Torlcvic h fc Melnicd 
Il98ll : iTenorio-Tagle. Mu fioz-Tufion fc Cox 1993I ). and pre- 
dicts that the lin e profile should be approximated by 
fviriai = \/ GM/R. This explanation has been successfully 
used to explain and calibrate th e R oc a and L oc a rela- 
tions found for giant H ll regions (|Terlevich fc Melnicklll98ll : 
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iMelnick. Terlevich fc MolesI llQSSL and explain the super- 
sonic widths foun d in NGC 604 (jMufioz-Tufion et"aLlll996l : 
lYang et al.lll996l ). Although NGC 1569 cannot be likened 
to a single giant Hll region, we can make a crude esti- 
mate of the expected virial broadening by taJting the mass 
to be the total H i mass of the galaxy (Mhi ~ 10* Mq; 
iMiihle et all |2005|) a nd an isophotal radius of ~4 kpc 
l|de Vaucouleurs et al.lll99ll : lMiihle et al.lbOOSi ') . This results 
in a broadening of cTviriai « 11 kms~^, in surprising agree- 
ment with our measured FWHMq. 

Thus, in conclusion we find that CI represents the tur- 
bulent ISM, where its motions result from a convolution of 
both the underlying gravitationally induced virial motions 
inherent in the gas and the general stirring effects of the 
starburst. It is also possible that an additional minor con- 
tribution to its width may originate from unresolved shell 
components along the line-of-sight. 



6.2.2 The broad component, C2 

We now turn our attention to the intriguing broad compo- 
nent (C2). In the following discussion, we will attempt to 
understand what the conditions of the emitting gas are, and 
what mechanism(s) are acting to give it such a broad width. 
Firstly, we can use Figs [13] and [14] to identify a number 
of trends to help constrain the physical properties of this 
component. Fig. 1131 shows that the group of C2 points can 
be split into two further sub-groups. The narrower group 
(FWHMsilSO kms~^ labelled C2-1 on the plot) is associ- 
ated with the regions of higher C2 flux (compare the central 
panels of Fig. [6] with Fig. [7]), whereas the close grouping of 
high-FWHM C2 measurements («280 kms"^ labelled C2- 
2) correspond (to a high degree of correlation) with the most 
intense CI emission (i.e. the two bright knots; see Fig. [9] 
right panel). Disregarding the red tail of the CI velocity dis- 
tribution discussed earlier, there is an average red ward off- 
set of only ~10 kms~^ between CI and C2, implying that 
C2 does not originate from gas with large-scale bulk mo- 
tions. Fig.[l4]clearly demonstrates that the broadest C2 lines 
are a lso the faintest (a common result; Im ufioz-Tufion et al.l 
1 19961 ). where the average ratio of flux(C2)/flux(Cl) over the 
FoV is 0.07. We can now assert that any theory explaining 
the origin of C2 must describe both how it exists over the 
entire FoV, and why its width is at once so highly correlated 
with the position of the bright knots seen in the narrow com- 
ponent, and so anti-correlated with its own flux. 

Mechanisms proposed in the literature that can account 
for line widths of >100 kms~^ in Hll regions are quite var- 
ied, particularly from studies of more distant systems. More- 
over, there are few studies that compare in terms of both the 
spatial-scale of the overall region (i.e. 30 Dor, for example, is 
clearly much smaller than the NGC 1569 starburst) and the 
spatial-resolution of the observations. This makes it harder 
to evaluate which mechanisms may be at work in NGC 1569. 

Unresolved expanding shells have again been cited as 
a possible cause of the broad component in a few studies, 
but it is difficult to envisage how smooth Gaussian profiles 
with widths of up to 300 kms^^ could be fully accounted 
for in this manner without showing some s igns of velocity 
splitti ng at our level of spectral resolution. iHeckman et al.l 
j 19951 ) observed Hq emission from NGC 1569 using two per- 



pendicularly oriented long-sits intersecting at SSC A, and in 
positions where they were able to fit multiple Gaussians to 
the line profile, they found two distinct kinematic systems. 
Lines near to Days were found to have widths of FWHM ~ 
30-90 kms^^, whereas broader lines with widths of up to 
150 kms~^ were only identi fied at radial ve locities of over 
200 kms~^ relative to Ways - iHeckman et al.l proposed that 
this second, high-velocity system is the origin of the broad 
widths (30-50 A at full-width zero-intensity) found in the 
summed Ha line profile of the area surrounding SSC A, since 
when integrating a larger spatial area, the discrete kinematic 
systems became blended togetheiQ. The level of broaden- 
ing due to unresolved components is therefore dependent 
on both spatial- an d spectral-resolution. This is illustrated 
well by the study of lMelnick et al.l (|l999l ). described above. 
Although their high spatial- and spectral-resolution obser- 
vations allowed them to resolve individual line components 
(shells) on scales down to 0.1 pc, they were still able to iden- 
tify a clear, underlying, broad component persisting in all 
sight lines. 

Can the effects of high-energy photons and fast- 
flowing winds from the nearby clusters explain the 
origin of C2? Over the years, much work has been 
done examining these effects using complex hy- 
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mass between the different hot and cold phases of the ISM 
can occur through four main processes: photoevaporation; 
conductively-driven thermal evaporation; hydro dynamic 
ablat ion; and turbulent mixing (for a review see [Pittardl 
I2OO6I ). Photoevaporation takes place in the presence of a 
strong ionizing radiation field, and conductive evaporation 
occurs when energetic electrons from a hot, surrounding 
medium, deposit their energy in the cloud surface. When 
the hot, tenuous medium is fast-flowing, as in the case of 
stellar or cluster winds, shearing between the high velocity 
flow and the cloud surface can set up a turbulent mixing 
layer that can drive further evaporation. The impacting 
wind can also physically strip (ablate) material off the 
clu mp surface, and entra i n it in to the flow (|Scaldll987l ). 

iBegel man fc Fabian] (|l990l ) described a model for a tur- 
bulent mixing layer that forms between the hot, inter-cloud 
gas and the surface layers of cool or warm clouds in near- 
pressure equilibrium. The UV radiation emitted from the 
mixing layer photoionizes the cold gas in the cloud sur- 
face giving rise to optical emission lines (in a similar way 
to a shock precursor) with charac teristic velocities of a frac- 
tion o f the hot-phase sound speed. [Slavin. ShuU fc Begelman] 
1 19931 ) expanded on this idea by modelling the cloud inter- 
face in detail. Including the effects of nonequilibrium ion- 
ization and self-photoionization allowed them to predict a 
number of optical, IR and UV line ratios, which they found 
to be in good agreement with the observed values at high 
latitudes in our Galaxy (for which this model was originally 



* In our data, we associate C3 wi th part of the high- velocity sys- 
tem described by]j 



Heckman et al 



and where detected, wc resolve 



it at a sufficient level for it not to be confused with the broad C2. 
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developed). ISlavin et al] (|l993l ) show a schematic of their 
model in their figure 1: hot gas flows past a sheet of cold or 
warm gas inducing Kelvin-Helmholtz instabilities and the 
build-up of a turbulent mixing layer at the interface region. 
Thermal conduction between the layers results in a turbulent 
cascade down to the dissipative level, allowing the efflcient 
exchange of energy between the media. Their model predicts 
strong optical [O ill] and Ha emission with non-Gaussian 
(broad-winged) line profiles caused by the high levels of tur- 
bulence. Strong far-UV emission from C iv and O vi is also 
predicted. These predictions are qualitat ively in good agree- 
ment with what we observe, however, the lSlavin et al.l model 
predicts [Olll]/H/3 and [Sll]/Ha nebular line ratios firmly 
in the shocked regime, far away from our observed ratios. 
A decrease in the 'maximum starburst' line threshold (indi- 
cating gas with a shocked emission component) due to the 
low metallicity of NGC 1569 may account for part of this 
discrepancy. 

It is interesting to no te that observations of CO (e.g. 
Falgarone fc Phillipslll990D CH+ fCrane Hegyi fc Lamber j 
199ll ) and 2.1 urn (e.g. ICeballe et al.lll986h in close-bv 
Galactic molecular clouds have also revealed broad line pro- 
files similar in shape to what we observe in our optical 
recombinat ion and forbidden lin es. T o explain these ob - 
servations, iHartouist et al.l l|l992h and iDvson et all (119951 ) 
also developed turbulent boundary layer models, and found 
these broad molecular line widths could again be reproduced 
through the presence of shear-induced turbulence, resulting 
from shocks in the cloud interface due to the vi scous coupling 
between the fast wind and the clump gas. iMelnick et al.l 
l|l999h were the first to suggest a connection with turbulent 
layers and erosion of dense gas clumps in an extragalactic en- 
vironment when discussing the persistent underlying broad 
line component identified from their high-resultion observa- 
tions of 30 Dor. 

We therefore find the most likely explanation of the 
origin of C2 is emission from a turbulent mixing layer 
formed at the surface of dense clouds resulting from the 
viscous coupling between the cool clump material and the 
hot, fast, winds from the surrounding clusters. Material is 
also removed from this layer through thermal evaporation 
and/or mechanical ablation, resulting in the mass- loading 
of the fiow. We know from X-ray studies that the hot, dif- 
fuse medium in NGC 156 9 has a characteris tic temperature 
around 0.7 keV (~10^ K: lMartin et al.ll2002l '). implying that 
the sound speed in this hot phase is ~500 kms~^. The line- 
widths we measure are only 60 per cent this value, mean- 
ing that the turbulent motions are subsonic with respect to 
the hot gas. This may provide a further explanation for the 
lack of shock-excited line ratios in our data, since it is ex- 
pected that photoionization in these central regions would 
overwhelm the signature of shock excitation anyway. 

Although a shear-induced turbulent mixing layer can 
theoretically attain a quasi-steady state when the energy 
fiux i nto the layer is e qual to the cooling rate per unit 
area (|Slavin et al.lll993h . the additional destructive effects 
of evaporation and/or ablation mean that gas clumps will 
only have a limited lifetime. How long could we therefore 
expect the clouds to survive? Unfortunately our knowl- 
edge of cloud disintegration is not yet developed enough 
to answer this question in detail. Simple hydrodynamical 
simulations of ablation predict that clouds enveloped in a 



supersonic (relative to the cloud temperature) wind have 
survival timescales typicafl y less than 1 Myr (jKlein et al.l 
ll994l : lMarcolini et aLlfioOSl ). Clearly this cannot be the case 
since warm ionized gas is observed in superwinds out to 
heights of many kpc (e.g. Shopbell fc Bland-Hawthornlll998l : 
lMartin|[l99i : iDevine fc Ballvill999l '). Further detailed mod- 
elling work is required. 

Whether or not these mechanisms are all acting in our 
particular FoV, we conclude that the energetics are domi- 
nated by the interface layer at the surface of the gas clumps, 
and that the principle infiuence on the gas state is external 
to the clouds, rather than embedded within the cloud, as 
would be the case for isolated giant Hll regions. In other 
words, the properties of the Ha knots in our FoV are tied 
to the intrinsic state of the gas rather than the available en- 
ergy supplies (the star clusters). This important difference 
may help to explain why the turbulent widths decrease away 
from the cloud surfaces. Cloud material remains intact and 
cool only for a limited length of time after it leaves the sur- 
face of the gas clump. Dense gas found near the surface of 
the cloud would be able to emit in Ha at higher velocities 
as it becomes entrained into the high-velocity flow than less 
dense gas, resulting in a broader line-width at the location 
of the cloud. 



6.2.3 The motion of C2 

The average width of C2 is much higher than the average 
radial velocity offset between CI and C2 (obvious from com- 
paring Figs [T3] and I14|l , meaning that in our FoV, turbu- 
lent m otions dominate over bulk flows. iTenorio-Tagle et al.l 
(|2006l ) develop 2D hydrodynamical models of the evolution 
of an Hll region assuming that the medium surrounding 
the ionizing star cluster is composed of a large collection of 
clouds (i.e. a clumpy medium). Although their model is not 
speciflcally applicable to our case since the ionizing source is 
located within the cloud stratum, it illustrates a number of 
interesting results regarding the flow of material within the 
outer bounding shock front. In their simulations, the lead- 
ing shock filters slowly through the gas clouds dissipating 
them one- by-one, leaving a highly turbulent, mass-loaded 
medium in its wake. The material interior to the leading 
shock moves in all directions resulting from the obstruction 
by the gas clumps, and possesses only a small net velocity 
in the outward direction. This is a very similar situation to 
what we observe, and shows that our FoV must be sampling 
gas well within the bounding shock (or sonic point) of the 
galactic wind, since as the flow becomes supersonic, turbu- 
lence is damped leaving only the bulk velocity component. 



6.2.i Final remarks 

We find small-scale variations in the velocity fields of both 
CI and C2 components, suggesting local energy inputs. 
However, these cannot (yet) be tied unambiguously to the 
effects of any single energy source (e.g. cluster 10) . The dis- 
tribution of non-photoionized emission covering the whole 
of the south-eastern portion of the FoV (see Fig. Ilip would 
suggest that high-velocity gas is impacting from the south- 
east (i.e. the direction of SSC A), but without information 
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on the three-dimentional distribution of the gas and velocity 
field, we cannot be sure. 

Since the constraints to our model of the origin of C2 
are only qualitative, further observations and modelling are 
clearly required. We predict that if the energetics are dom- 
inated by the interface region where the turbulence is pro- 
posed to originate, observations of emission originating from 
deep within the cloud would not show this line broadening. 
This would also apply to regions shielded from the wind 
impact, where only the ambient IS M turbulence would be 
evident. Furthermore, the models of iMarcolini et al] (|2005l ) 
predict strong soft X-ray emission from the bow-shock (re- 
verse shock) produced as the wind hits the cloud. Whether 
the X-rays are primarily emitted from wind material or 
evaporated cloud material depends on the metallicity of the 
wind and the strength of the thermal conduction within the 
cloud. UV emission (or absorption) from C iv and O vi is 
also predicted to arise through collisional processes within 
the cloud m ateri al at the wind- cloud interface by both 
ISlavin et al.l and IMarcolini et al.l Currently, O vi-emitting 
gas is the highest temperature phase of an outflow that can 
be probed for which r eliable velocities can also be measured 
l|Heckman et aLlboOll ) , making this transition a very impor- 
tant tool in the investigation of galactic winds. 

In a companion paper (Westmoquette et al. 2007a; Pa- 
per II) , we present the analysis of the other three IFU fields 
covering the central regions of NGC 1569. Following on from 
our analysis presented here, we first concentrate on the de- 
tails of each field, then combine all the data in an attempt to 
constrain the general state of the ISM in the whole region. 



7 SUMMARY 

The unique capabilities of the CMOS IFU, coupled with 
HST imaging, have allowed us to map the detailed ionized 
gas properties of a 5 x 3.5 arcsecs region near the nucleus of 
NGC 1569, and to relate the presence of a young star cluster 
located in the centre of the field-of-view to the surrounding 
gas. This field is part of a large, bright Hll region complex 
(Hll region 2: IWalleilll99ll ) that dominates the emission in- 
tensity, and exhibits evidence for recent, on-going and cm- 
bedded star- and cluster-formation, presumab ly fuelled by 
the l arge near-by repository of molecular gas (jTavlor et al.l 
1 19991 ). 

Summing the spectra from an aperture centred on the 
position of the cluster has allowed us to measure and derive 
a number of fundamental characteristics of this object. We 
find: 

• Cluster 10 is composed of two sub-clusters, which we 
have denoted lOA and lOB, with sizes -Roff.A = 0.88 and 
-Reff.B = 0.60 pc and a projected separation of 3.7 pc. Be- 
cause of their proximity, they cannot be resolved in our IFU 
data. 

• Using model colour-colour diagrams, we were able to 
estimate their ages to be 5-7 Myr for 10 A and $;5 Myr 
for lOB, with reddenings consistent with the Galactic fore- 
ground level (confirmed by measurements of the Ha/H/J ra- 
tio from our IFU spectrum). 

• We derive a photometric mass of 7 ± 5 X 10^ Mq from 
the combined light of both sub-clusters. This result means 



that cluster 10 is not a super star cluster in the generally 
accepted definition. 

• The presence of the red WR-bump (resulting from 
broad C iv emission) confirms the young age of at least one of 
the sub-clusters. From measuring the total luminosity of the 
feature, we estimate an equivalent population of 10 ±5 WC4 
stars in cluster 10. This is consistent with the narrow-band 
Hell observations of iBuckalew et al] (|2000l ) who find emis- 
sion associated with lOB equating to a luminosity equivalent 
of 3 WNL stars. 

The spatially resolved nature of our observations have al- 
lowed us to make a detailed analysis of the distribution of 
gas within the FoV. The most important findings are: 

• We find bright emission over the whole field well char- 
acterised by a Gaussian function with mildly supersonic 
line widths (CI), and a distinct lower limit of FWHMq ~ 
35 kms~^ (equating to a = 12 kms^^ after correction for 
thermal broadening). 

• We also measure a fainter line component (C2) with 
line widths of between --^100-300 kms~^, also found over 
the whole field. 

• The width of the broad underlying emission is both 
highly correlated with the position of the bright narrow 
component emission (i.e. bright CI oc broad C2), and anti- 
correlated with its own fiux (i.e. bright C2 cx narrow C2). 

• We conclude that CI represents the turbulent ISM, 
where its motions result from a convolution of gravitational 
virial motions and the general stirring effects of the star- 
burst (giving rise to the clear lower-limit to the line width, 
FWHMo). However, it is possible that an additional minor 
contribution to its width may originate from unresolved shell 
components along the line-of-sight explaining the local vari- 
ations in CI width across the FoV. 

• We find the most convincing mechanism to explain the 
broad component is emission from a turbulent mixing layer 
on the surface of the cool gas clumps in the ISM of NGC 
1569. This layer results from the coupling between the clump 
material and the hot, fast winds from the surrounding clus- 
ters, and gives rise, through thermal evaporation and/or ab- 
lation, to a tenuous, highly turbulent, diffuse gas that per- 
vades the local environment that will eventually mass-load 
the outflow. 

• The distribution of non-photoionized (shocked) gas as 
determined from nebular diagnostic line ratios supports ev- 
idence that the dominant source of wind energy /material is 
from SSC A, consistent with morphological indicators such 
as swept-back cones around dense Ha emitting knots that 
point towards this cluster. This does not mean, however, 
that the signiflcant proportion of ionizing flux in this region 
does not come from the young cluster 10. 

• The average radial velocity offset between CI and C2 
is very small compared to the average line- width of C2, in- 
dicating that in our FoV, turbulent motions dominate over 
bulk flows, and that we are sampling well within the 'energy 
injection zone' of the outflow. 

Our results reflect the presence of a complex ISM with 
a wide range in physical conditions. Our suggestion that gas 
is being mixed from the relatively cool photoionized medium 
into the hot phase carries with it a number of astrophysical 
implications that could generally apply to galactic winds. 
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Among these would be mass loading of the hot ISM occurs 
close to or within the star forming disk. This will enhance 
the hot ISM cooling rate and increase the density of material 
injected into the base of the thermal wind, thus influencing 
the properties of the wind. A full discussion of the ionized 
gas state in NGC 1569 will be presented in a companion 
paper (Westmoquette et al. 2007a; Paper II), where we put 
our findings for this IFU position in context with the whole 
central region using observations made with three other IFU 
pointings. In a forthcoming contribution (Westmoquette et 
al. in prep.; Paper III) we will examine the properties of 
the galactic outflow as a whole using deep Ho? imaging and 
IFU observations and attempt to create a unified physical 
picture of the outfiow state. 
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